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I Abstract. In this paper, we first prove the unique global strong solution with 

vacuum to the two dimensional nonhomogeneous incompressible MHD system, 
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as long as the initial data satisfies some compatibility condition. As a corollary, 
the global existence of strong solution with vacuum to the 2D nonhomogeneous 
incompressible Navier-Stokes equations is also established. Our main result 
. improves all the previous results where the initial density need to be strictly 

positive. The key idea is to use some critical Sobolev inequality of logarithmic 
type, which is originally due to Brezis-Wainger [7]. 
'■""J ' Keywords: nonhomogeneous incompressible fluid, strong solution, vacuum. 
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1. Introduction 

Magnetic fields influence many fluids. Magnetohydrodynamics(MHD) is con- 
cerned with the interaction between fluid flow and magnetic field. The governing 
equations of nonhomogeneous MHD can be stated as follows [13] , 
(1.1) 

' p t + div(pw) = 0, in Q. x [0, T), 

(pu) t + div(pu ® u) - div(2p(p)d) - (B ■ V)B + VP = 0, in Q x [0, T), 

B t - XAB - curl(u x B) = 0, in Q x [0, T), 

k divw = 0, dwB = 0, in Q x [0, T). 

Here p and u are the density and velocity field of fluid respectively. P is the 
pressure. B is the magnetic field. p(p) > denotes the viscosity of fluid, which 
we assume in this paper is a positive constant. A > is also a constant, which 
describes the relative strengths of advection and diffusion of B. For simplicity of 
writing, let p, — A = 1, d = | (Vit + (Vu)') is the deformation tensor. 

In this paper, we focus on the system ( II .ip with the initial-boundary conditions 

(1.2) u = 0, B-n = 0, curLB = onffix [0,T), 

(1.3) (p,u,B)\ t=0 = (p ,u ,B ) inn. 
Here Q is a bounded smooth domain in M 2 . 

If there is no magnetic field, i.e., B = 0, MHD system turns to be nonhomoge- 
neous Navier-Stokes system. In fact, due to the similarity of the second equation 
and the third equation in (II. ip . the study for MHD system has been along with 
that for Navier-Stokes one. Let's recall some known results for 3D nonhomo- 
geneous Navier-Stokes equations. When the initial density po is bounded away 
from 0, the global existence of weak solutions was established by Kazhikov |21] . 
see also [3]. Moreover, Antontsev-Kazhikov-Monakhov [S] gave the first result on 
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local existence and uniqueness of strong solutions. For the two-dimensional case, 
they even proved that the strong solution is global. But the global existence of 
strong or smooth solutions in 3D is still an open problem. For more results in 
this direction, see P2H2UI2E] an d references therein. 

If the initial density po allows vacuum, the problem becomes more complicated. 
Simon [29] proved the global existence of weak solutions, see also [26J. Choe- 
Kim [12J constructed a local strong solution under some compatibility conditions 
on the initial data. More precisely, they proved that if (po, ^o) satisfy 

(1.4) < po G L*(fi) nH 2 (n), n GDj(fi)nZ) 2 (fi), 
and the compatibility conditions 

i 

(1.5) divu = 0, — pAu + VPo = Pl9-, hi fl, 

with some (Po, g) belonging to D 1 (Q) x L 2 (Q), then there exists a positive time T 
and a unique strong solution (p,u) G C([0, T); H 2 (Q)) x C([0, T); D£(fi) nD 2 (Q)) 
to the nonhomogeneous Navier-Stokes equations, where Dq(Q) and D 2 (Q) denote 
the usual homogeneous Sobolev spaces. Recall that -Dq(R 3 ) = {u G L 6 (IR 3 ) : 
Vu G L 2 {R 3 )} and Dj(fi) = flj(fi) if Q CC M 3 . 

After the local existence of strong solution, one question came out naturally, 
which is whether the solution blows up in finite time. Suppose the finite blow-up 
time T* exists, [22J proved the Serrin type criterion, which says that 

f T * 2 n 

(1.6) J \W(t)\\ s L ^dt = oo, for any (r, s) with — I — = 1, n < r < oo, 

where n is the dimension of the domain and L r w is the weak U space. (The 
proof was given in [22] only for 3D case, but almost the same proof works for 2D 
case.) In particular, for the 2D case, it follows from the energy inequality the 
solution satisfies that sup 0<T<T »(|| v /pM|| L oo( 0T;i 2) + || Vu|| i 2( 0T;i 2)) is bounded, 
which implies that u G L 4 (0,T*;L 4 ) if p is bounded away from 0. Hence the 
criterion ( 11. 6ft in fact implies global existence of strong solution provided that po 
is bounded away from 0. However, if the density is allowed to vanish, whether 
the strong solution exists globally remains unknown. This is the main problem 
we shall address in this paper. 

Let's go back to the MHD system (II. ip . As said before, the research for MHD 
goes along with that for Navier-Stokes equations. The results are similar. When 
p is a constant, which means the fluid is homogeneous, the MHD system has 
been extensively studied. Duraut-Lions [17] constructed a class of weak solutions 
with finite energy and a class of local strong solutions. In particular, the 2D 
local strong solution has been proved to be global and unique. While for the 
three-dimensional case, different Serrin type criteria similar to ( II. 6p were given 
in [8l[T9j[20l|30] . As for the 3D Navier-Stokes equations, whether the local strong 
solution is global is still open. 

When the fluid is nonhomogeneous, Gerbeau-Le Bris [16], Desjardins-Le Bris 
[T4"] studied the global existence of weak solutions of finite energy in the whole 
space or in the torus. Global existence of strong solutions with small initial 
data in some Besov spaces was considered by Abidi-Paicu [TJ. Moreover, [1] 
allowed variable viscosity and conductivity coefficients but required an essential 
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assumption that there is no vacuum (more precisely, the initial data are closed to 
a constant state). Chen- Tan- Wang [TU] extended the local existence in presence 
of vacuum. In conclusion, if the initial data satisfies that 

(1.7) 0< Po eH 2 , {u ,B )eH 2 , 
and the compatibility conditions 

Mo = 0, Bq ■ ft — 0, curl£>o = 0, on dQ, 

(1.8) j 

divuo = div^o = 0, -Au + VP - (B ■ V)B = p*g, in ft, 

with some {Po,g) G H l x L 2 , then there exist a positive time T and a unique 
strong solution (p,u,B) to the problem fll.lD - fll.3p . such that 

peC([0,T];H 2 ), (u,B)eC([0,T];H 2 ), 

(1.9) p G C([0, T};H l ) n L 2 (0, T; # 2 ), K 5*) G L 2 (0, T; tf 1 ), 

and (p t ,^pu t ,B t ) GL°°(0,T;L 2 ). 

For all the techniques, refer to [TT] . 

It comes to the question whether the local strong solution blows up. After the 
proof of [22J for nonhomogeneous Navier-Stokes equations, one can get the same 
criterion (jl.6p for nonhomogeneous MHD, see also |31] . In particular, for the 
2D case, it says that ||m||l2 L oo becomes unbounded once the local strong solution 
blows up. On the other hand, the energy inequality tells us || Vm|| L 2 L 2 is uniformly 
bounded, which only imply that ||«||£2( BAf0 _) is uniformly bounded. Therefore, 
in view of the blowup criterion ( II. 6p . it's not enough to extend the local strong 
solution to global one. To improve the regularity of the velocity, we choose to 
apply a critical Sobolev inequality of logarithmic type, which is originally due to 
Brezis-Gallouet [6] and Brezis-Wainger [7j. In this paper, we use some extension, 
which was proved by Ozawa [27]. For a new proof, see |23j. The inequality is 
stated as follow, 

Lemma 1.1. Assume f G H 1 ^ 2 ) n W^iR 2 ), with some q > 2. Then it holds 
that 

(1.10) imU°o (R2) < C (l + ||V/|| La(Ra) (ln + WfWw^Q*))*) , 

with some constant C depending only on q. 

The same proof with some proper extension theorem(see [2]), in fact gives the 
following modified inequality, which involves the integral with respect to time. 
For completeness, we will give the proof in Section 2. 

Lemma 1.2. Assume Q is a bounded smooth domain in M. 2 and f G L 2 (s, t; iif 1 (fi))n 
L 2 (s,t; W 1,q (Q)), with some q > 2 and < s < t < oo. Then it holds that 

(1-11) ||/||£3( a ,ti£»(n)) <c(l + ||/||ia( s ,t;fli(n)) (hi + WfW^is^w^m)^) » 
with some constant C depending only on q and Q, and independent of s,t. 

The application of ( 11. lip is the key idea of this paper. Due to this, we can close 
the estimates for \\(u, B)\\L<x>gi. The higher order estimates are in the same spirit 
of [22]. For more details, see Section 3. Finally, we get the result about global 
existence of strong solution. 
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Theorem 1.3. Assume that the initial data (pq,Uq, B ) satisfies (ll.7j) and the 

compatibility conditions (jl.8p . Then there exists a global strong solution (p,u,B) 



of the MED system (Oil - (Q)l . wito 

p G C([0, oo); if 2 ), (u, B) G C([0, oo); tf 2 ), 
(1.12) P G C([0, oo); H 1 ) n L 2 OC (0, oo; if 2 ), (w<, 5*) G L 2 OC (0, oo; tf 1 ), 

and (p t , ^ut, Bt) G L~ (0, oo; L 2 ). 



Some remarks are given about this theorem. 

Remark 1.1. The local existence of unique strong solution with vacuum to the 
system II 1.1}) in a two-dimensional bounded domain can be established in the same 
manner as [IB] and flU/ . Through this paper, we will concentrate on establishing 
global estimates for the density, velocity and magnetic field. 

Remark 1.2. If we consider the most special case, where p is a constant(the fluid 
is homogeneous) and B = 0(no magnetic field), then the system (11.11) becomes 
the classical Navier-Stokes system. The global existence of strong solution has 
been proved by Leray 125^ . More generally, if we consider the case that only p is 
a constant, the system (II. ip becomes the classical homogeneous MHD system. As 
said before, the corresponding result has been derived by Duraut-Lions [T 7 ^ . 

If B = 0, Theorem 11.31 in fact gives a positive answer to the global existence of 
strong solutions with vacuum of the 2D nonhomogeneous Navier-Stokes system. 
It covers the corresponding result in [5], where the density is strictly positive. 

Corollary 1.4. Assume that the initial data (po,wo) satisfies (11.71) and the com- 
patibility conditions ( 11.51) . Then there exists a global strong solution (p,u) of the 
Navier-Stokes equations, with 



We conclude this section with some notations and lemmas. Z/(f2), W k,r (fl), 
(1 < r < oo), are the standard Sobolev spaces, and we use U = L r (Q), W k,r = 
W k ' r (tt). Especially, when r = 2, denote H k = W k ' 2 . For simplicity, let 



Some more lemmas will be used during the proof of Theorem 11.31 One is 
following from the regularity theory for Stokes equations. For its proof, refer 
to PSJ. 

Lemma 1.5. Assume that (u, P) G Hi x H 1 is a weak solution of the stationary 
Stokes equations, 



(1.13) 



p G C([0, oo); H 2 ), u G C([0, oo); H 2 ), 

P G C([0, oo); H l ) H L 2 OC (0, oo; H 2 ), u t G L 2 OC (0, oo; H 1 ) 

and (p t , y/pu t ) G L% c (0, oo; L 2 ). 




(1.14) 




Au + VP = F, inVL 
r u = 0, in Q, 
= 0, on dQ, 



;i.i7) 
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and FaL q ,l<q<oo. Then it holds that 

(1.15) IMIw 2 '8 < CH-^llw + C||u||^x, 

with some constant C depending on Q and q. Moreover, if F 6 H 1 , then 

(1.16) \\u\\ H s < C\\F\\ H i + C\\u\\ H i, 
with some constant C depending only on Q. 

The other lemma is responsible for the estimates for B and follows from the 
classical regularity theory for elliptic equations. For its proof, refer to [3]. 

Lemma 1.6. Assume that B e H l is a weak solution of the Poisson equations 

AB = G, in fi, 
B ■ n — 0, curlf? = 0, on d£l, 

and G^L q ,l<q<oo. Then it holds that 

(1.18) \\B\\ W 2 :q < C\\G\\ Lq + C\\B\\ Hh 

with some constant C depending on Q and q. Moreover, if G G H l , then 

(1.19) \\B\\ H z < C\\G\\ H i + C\\B\\ H i, 
with some constant C depending only on Q. 

2. Proof of Lemma 1.2 

This section is dedicated to the proof of Lemma 1.2. First we will prove the 
inequality ( II. lip for the whole space case, which is 

(2.20) \\f\\L 2 (s,t;L°°(M. 2 )) < C (l + ||/||x2(s,t;H 1 (K 2 )) (hl + ||/|U 2 («,t;TVl.9(R 2 ))) 

The proof follows exactly that in [23J and lies mainly on the Littlewood-Paley 
decomposition. So we introduce here some new notations associated with the 
decomposition. Define C to be the ring 

C = {^K 2 : f<KI<f 

and define T> to be the ball 

V = ^eR 2 : lei < \ 

Let x an d <p be two smooth nonnegative radial functions supported respectively 
in T> and C, such that 

X(0+X>( 2 "*0 = 1 for^GM 2 , and £>(2-*0 = 1 for ^ G M 2 \ {0}. 

<?GN q& 

Denote the Fourier transform on M 2 by F and denote 

h = F~ l tp, h = F^x- 
The frequency localization operator is defined by 

AJ = F- 1 [p{2- q i)FU)] = 2 29 / h(r y )f(x - y)dy, 
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and 

SJ = J"' 1 [x{2-^)Hf)\ = 2 29 f h(2«y)f(x - y)dy. 
Now it's ready to prove (I2.20p . 
Proof. Decompose / into three parts such as 

f(x, t) = S-M-ifix, r)+J2 r) + A if( x > T ) 

(2.21) |j|<jv j>N 

= ft(x, r) + f 2 (x, r) + f 3 (x, t). 
By Bernstein's inequality (see [5]), 

(2.22) H/ilU^-) < C- 2N ^\\f\\ LHs ,, Lqh q G [l,oo). 
Similarly, 



11/2 1| L 2 (s,t;L°°) < ^2 \\^if\\L 2 (s,fr,L°°) 
\j\<N 



(2.23) 



and 



<cm (||v(A,/)||i 2(Sit;i2) 

<CiV5||v/|| i2(Sit . i2) , 



||/3||L 2 ( s ,t;L^) < II ^jf\\L 2 (s,t;L°°) 

j>N 

(2-24) <L7^2 2 ^i/ 2 )|| V/ || L2(M;L9) 

j>N 

If we set K = min(2/g, 2(1/2 — 1/q)), then 
(2.25) ||/ ||z W ~) < C {2- k7V ||/|| l2(Si ^ 1 „ ) + Nh \\ V f\\ LHs>t;L2) } . 



Choose N 



log 2 > 



11/11 

l|V/|| i a ( . lt . x ,2 ) J 



+ 1, hence we derive that 



(2.26) \\f\\ L2(s ,, L ^ < C7||V/||^(^) ( 1 +(^ + l^n a ' t ^ 1 ' q) ) 

\ V l|V/||L2( S ,t;L2)/ 

which implies fl2T20|) . □ 
Combining the extension theorem(see |2]) and (12.201) . we prove Lemma 1.2. 

3. Proof of Theorem 1.3 

This section is dedicated to the proof of Theorem 1.3. Define the quantity 
$(T) as follow, 

$(t)= sup (1^)11^ + \\um^ + \\Bm^) + \\vp^\\l-(o,T^) 

0<t<T 

(3.27) pT pT 

+ {\Ht)\\ 2 HS + \\B(t)\\ 2 H s)dt+ {\\u t \\ 2 Hl + \\B t \\ 2 Hl ) dt. 
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Suppose the local strong solution blows up at T* < oo, we will prove that in 
fact there exists a generic constant M < oo depending only the initial data and 
T* such that 

(3.28) sup $(T) < M. 

0<T<T* 

Having (13. 28j) at hand, it is easy to show without many difficulties that we can 
extend the strong solution beyond T*, which gives a contradiction. Hence the 
local strong solution does not blow up in finite time. Also, the uniqueness of 
strong solutions is a standard procedure. 

Through out this section, C denote a generic constant only depending on the 
initial data and T*. The proof is divided into five steps, due to different level 
estimates. 

Before proceeding, we write another equivalent form of (II. ip for convenience, 
which is 

( p t + u ■ Vp = 0, 
pu t -Au + (pu ■ V)u - (P ■ V)P + VP = 0, 
B t - AP + (u ■ V)P - (B ■ V)« = 0, 
divw = 0, divP = 0. 



(3.29) 



Now we start the proof of Theorem 11.31 

Step I L°° bound for p. The equation ( I3.29p i for density is a transport 
equation, then for every < t < T*, 

(3.30) IIp(*)IU» = IIa)IU- 

Step II Basic energy estimate 

Proposition 3.1 (Energy inequality). There exists a constant M depending only 
on \\y/p^u \\ L 2 and ||P ||l 2 ; such that for every < T < T* , 

(3.31) ||v / P«ll!-(0,T;L 2 ) + ll^lli- ( o,T;L 2 )+ / II Vufpdt + f \\V B\\\,dt < M . 

Jo Jo 

Proof. The proof is standard. Multiplying (I3.29l) 2 and (I3.29[) 3 by u and B respec- 
tively, then adding the two resulting equations together, integrating over Q, one 
can get that 

(3.32) ^/ /0 l M l 2rfx + ^/l jB ' 2rfx + J\^ u \ 2dx + f\VB\ 2 dx = 0, 

where integration by parts was applied. It implies that the inequality ( 13 . 3 1 [) holds 
and consequently completes the proof. □ 

Step III Estimates for \\(y/pu t , Pt) ||l 2 (o,t ; l 2 ) and ||(Vw, VP)||l°o( 0iT;L 2). 

This is a crucial step during the proof. Higher order estimates of the den- 
sity, velocity and magnetic field can be done in a standard way provided that 
\\(u, B)\\ H i is uniformly bounded with respect to time. To prove that, we will 
make use of some extension of critical Sobolev inequality of logarithmic type, as 
indicated by Lemma [1.21 
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Proposition 3.2. Under the assumptions in Theorem \1.3[ it holds that 

(3.33) sup B(T))\\ 2 H1 + f \\{^pu u B t )tf»dt\ < oo. 

0<T<T* L JO J 

Proof. Multiplying the equation ( I3.29p 2 by u t and integrating over Q lead to 

(3.34) --^ J \Vu\ 2 dx + J p\u t \ 2 dx = — J \pu • Vw) • u t dx + J {B -V)B ■ u t dx. 
By Holder's inequality and Young inequality, 



(3.35) 



(pu ■ V)w • u t dx 



< C\\yfpU t \\ L 2 ■ \\u\\ L ™ ■ ||VW|| L 2 



< ^Iv^ll^ + ^NIi-IIVu^- 



Applying integration by parts with the conditions that divl? = in Q and B-n = 
on dQ, then 



(3.36) 



d 
"dt 



[B ■ V)B ■ u t dx 

(B ■ V)B -udx- J (B t ■ V)B -udx- J {B ■ V)B t ■ udx 

-± j(B. V)w • Bdx + J(B t - V)m • Bdx + J (B ■ V)u • B t dx 

d f 1 
<-- (B-V)u-Bdx + C\\B\\ 2 L ^\\Vu\\ 2 L , + -\\B t \\ 2 L 2. 

Hence, combining ( 13.34[) -( l3736|) . we get that 

I 1 d f d f 

2^ Ut ^ 2 + 2lt J lVul2dX + Jt j( B - V ) u - Bdx 



(3.37) 

Similarly, multiplying the equation ( I3.29p 3 by B t and integrating over Q lead 



to 



1 d 
2~dt 



\VB\ 2 dx 



\B t \ 2 dx 



(3.38) 



- / (u ■ VB) ■ B t dx + / (B ■ V)w ■ B t dx 



<2ll^ll^ + ^lkll^livs|| 2 i2 + c||v w ||| 2 ||s||! 



which implies that 

(3.39) — J \WB\ 2 dx+ \\B t \\ 2 L 2 < C||u||i,»||VS||ia + C\\B\\ 2 LOO \\\7u\\ 2 L2 . 

The term J (B ■ V)u ■ Bdx on the left hand of (13.371) can not be determined 
positive or negative, so we choose some appropriate positive terms to control it. 
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Note that it follows from Gagliardo-Nirenberg inequality that 



(B ■ V)u ■ Bdx 



<\\B\ 



L 4 



(3.40) 



<C\\B\ 



l 2 



\B\ 



L 2 



H 1 



iVd 



L 2 



<4l|Vtt||i, + C 1 ||S||i a (||S 



I 2 

Il-> 



\VB\ 



L 2 i 



Next, we multiply (13.391) by 2C\M + 2, where C\ and M are constants appearing 
in (13.401) and (I3.3ip . add it to (I3.37P and integrate with respect to time, then for 
every < s < T < T*, 

(3 ' 41) 

J \Vu(T)\ 2 dx + J \VB(T)\ 2 dx + J WVputWhdr + J \\B t \\ 2 L2 dr 

expjc^ ( || u || \ x + \\B\\ 2 Loo )dr^ + C. 



<C 



\Vu(s)rdx+ / VB(s)\ 2 dx 



Denote 



(3.42) *(f) = e+ sup (|k(r)||| 1 + || J B(r)|| 2 Hl )+ [ + \W\h) *r> 

0<T<t Jo 

then (EUD and fl33Tj) give that for every < s < T < T*, 



T 



*(T) < Ctf (s) exptC / (IMIioc + \\B\\ 2 Lao )dr 



(3.43) 



To get a proper estimate for ||m|| L 2 L oo and H-BH^^oo, we get help from Lemma 

o 

(3.44) 

IMIl 2 (s,T;.L°°) + II-^IIl 2 (s,T;L°°) 

<C|l + (||u||£ a( , ir . jr i ) + ll-Bll^^.jj!)) (ln + \\u\\ L 2 {SjT . w i A) +ln + H^H^^M)) } 
Applying Lemma [L5l to the equation (I3.29L yields 

(3.45) \\u\\ w x,4 < C\\u\\ m + C\\pu t \\j + C\\(pu ■ V)u - (B ■ V)B\\ L $, 

which implies 
(3.46) 

||m||l2( SjT;W i,4) <C||u||£,2( SiT . ff l) + C\\y/pUt\\ L 2 (s,T;L 2 ) 

+ ClklU^s.T;^ 1 )!! V-u|| L oo( S)T . L 2) + C||5|| L 2( S)T .^i)|| V-B||l°°(s,T;L 2 )- 



Similarly, applying Lemma [TTBl to the equation (I3.29IL to obtain 
(3.47) 

II-SIU 2 ^;^ 1 . 4 ) <C\\B\\ L 2^ sT . Hl ^ + C\\B t \\l 2 (s,T;L 2 ) 

+ C\\u\\ L 2( SjT . H i)\\V '-B||l°°(s,T;L 2 ) + C||-B|| L 2( S)T;H i)|| Vm||loo( SiT;L 2). 



Note that the constant C in (I3.46j) and (I3.47P does not depend on u, B, s or 
T. It only depends on the domain Q. Taking the energy inequality (13.311) into 
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consideration, then for every < s < T < T* 
(3.48) 



I M IIl 2 (s,T;L°°) + II^IIl 2 (s,T;L°°) 



<C 2 {l + (\\u\\l Hs ^ Hl) + \\B\\l HstT . H , } )\n(C(M,T*MT))) , 

where C 2 is constant which only depends on Q, and C(M,T*) is a constant 
depending on M in (I3.3ip and T*. 

Substituting (I3.48j) into (I3.43p . it arrives at 

(3.49) #(T) < C#(s) [C(M, T*)\I>(r)] C ' 2 (" u "^ 2 ( s ' T ' Hl ) + "' B "^ 2 ( s ' T: - ff1 )) . 

Recall the energy estimate (I3.3ip . one can choose s close enough to T*, such that 

(3-50) \im t C 2 (IH^t.h^ + \\B\\ 2 L 2 {SjT . H r^ < -, 

then for every s < T < T*, we have 

(3.51) q{T) <C^{s) 2 -C{M,T*) 2 , 

which completes the proof of Proposition 13.21 □ 

Remark 3.1. Unfortunately, we can not get any explicit bound for \\(u, B)\\ H i 
in terms of the initial data, due to the technique used here. 

We have some more estimates as corollaries of Proposition 13.21 
Proposition 3.3. Assume that 

(3.52) sup \uu(T),B(T))f Hl + f ||(V^ t , B t )f L *dt\ < C 3 . 

0<T<T* L JO J 

Then there exists a constant C4 depending on C%, such that 

(3.53) sup {||w||l2(o,t ; h2) + ||-B|| L 2( 0jT;H 2)} < C 4 . 

0<T<T* 

Proof. The equation (I3.29l) 2 . together with Lemma [1. 5 [ gives us that 

\\u\\ H * < C\\u\\ H i + C\\pu t \\ L 2 + C\\(pu ■ V)u\\ L 2 + C\\(B ■ V)B\\ L 2 

< C\\U\\ H 1 + C\\y/pU t \\ L 2 + C||m|| L oo||Vm|| L 2 + C||S|| L cx,||V5|| L 2. 



Similarly, by Lemma [1.61 

(3.55) \\B\\ H 2 < C\\B\\ H i + C\\B t \\ L 2 + C\\u\\ L <» || VB\\ l * +C\\B\\ l «,\\Vu\\&. 
Combining the two inequalities (13.541) and (13.551) . we have 

\\U\\ H 2 + \\B\\ H 2 

<C\\y/pUt\\v + C\\B t \\ L 2 + C (\\u\\ L oo + \\B\\ L oo + 1) • (\\ U \\ H 1 + \\B\\ H l) 

(3 56) 

<C(\\u\\ H 2 + \\B\\ H 2) 1 ' 2 {\\u\\ L 2 + ||i?|| L 2) 1/2 ■ {\\u\\ m + \\B\\ m ) 
+ C{\\u\\ m + \\B\\ m ) + C\\^-pu t \\ L 2 + C\\B t \\ L 2. 
where Gagliardo-Nirenberg inequality was used. Hence, 
(3.57) \\u\\ H 2 + \\B\\ H 2 < CWyfputWv + C\\B t \\ L 2 + C(l + \\u\\ m + H^M 3 , 
which completes the proof for (I3.53p . □ 
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Proposition 3.4. Assume (13.521) holds, then there exists some constant C 5 de- 
pending on C3 such that 

(3.58) sup {||m||l4(o,t ; l°°) + ||-B|| L 4( 0jT . £OO )} < C 5 . 
o<t<t* 

Proof. By Gagliardo-Nirenberg inequality, 

(3.59) \\uWloo < C||u||^ 2 ■ \\u\\^l 
and 

(3.60) \\B\\ L oo < C\\B\\% 2 ■ \\B\\^ 2 , 

which together with ( I3.53|) completes the proof for (I3.58p . □ 

Step IV Estimates for \\(^/pu t , B t ) \\l°°{o,t;V) and ||(Vw t , VB t )\\ L ^T;L^) 
From now on, the estimates are standard, due to the proof in [22J. We write them 
down here for completeness. 

Proposition 3.5. Under the assumptions in Theorem \1.3l it holds that 

(3.61) sup \\\(^pu t (T), B t {T))\\ H i+ I \\{Vu u VB t )\\ 2 L2 dt) < 00. 
o<r<r* L Jo J 

Proof. Taking t-derivative of the equation (I3.29L . then one gets that 
pu tt + {pu ■ V)u t - Au t + VPt 

(3 62) 

1 ' ; = -Ptu t - {p t u ■ V> - {pu t ■ V> + (B t ■ V)B + (B ■ V)B t . 

Multiplying (I3.62j) by u t and integrating over Q, 

— — f p\u t \ 2 dx + ( \Wu t \ 2 dx = — [ p t \u t \ 2 dx — I ' (ptu ■ V)w • u t dx 

(3.63) 2dt J J J J 

- I (pu t ■ V)u ■ u t dx + (Bf V)B ■ u t dx + (B ■ V)B t ■ u t dx. 

We estimate the terms on the right hand one by one. Taking f ll.ip i into consid- 
eration, we get that 

— J p t \u t \ 2 dx = J div (pu)\u t \ 2 dx 

(3.64) = - / 2pu ■ Vu t ■ u t dx 



< l\\v^\\h + c\\Vm\\hh\\l~, 



and also for the second term, 

— J (p t u ■ V)w • u t dx 

(3.65) = - pu- V[(w • V)w • u t ]dx 



< J \ P u t \M^ dx + / \ P u t \M^Wdx + / ,|«|»|v«||v«,|<i« 
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Here by Gagliardo-Nirenberg inequality, 

\put\ \u\ | Vu\ 2 dx 

(3.66) < llv^lU 2 IMU°°l|Vu|||4 

< C\\^u t \\ L 2\\u\\ La o\\Vu\\ L 2\\Vu\\ H i 

< Ni-iiv^tllia + ciivuii^i 

By Young inequality, 

I |/out||«| 2 |V 2 «|da; 

(3.67) 



a 



H 2 ' 



<C\\^pu t \\ L 4u\\ 2 La o\\V 2 u\\ L 2 

< h\\U\V^t\\h + c\\u\\ 2 H3 . 

And similarly, 

J p\u\ 2 \Vu\\Vu t \dx 

(3.68) <C\\u\\ 2 Loa \\Vu\\ L 4Vu t \\ L , 

<gl|V« t ||i a + C7||tt||l.||Vri||i 9 . 

For the third term of the right hand of f 1 3 . 6 3 [) . by Poincare inequality and 
Gagliardo-Nirenberg inequality, 

— j {pu t ■ V)w • Utdx 

(3.69) <C\\^pu t \\ L 2\\Vu\\ L 4u t \\ L 4 

<C\\u\\ 2 H 4^u4 2 L . 2 + l\\Vu t r L2 . 

Since divB t = in Q and B t ■ n = on dQ, then 

J (B t ■ V)B ■ u t dx 

(3.70) = - J(B t ■ V)u t ■ Bdx 

<l\\Vu4h+C\\B\\U\B t \\ 2 L2 . 



And similarly, 
(3.71) 



J(B-V)B t - u t dx 



<l\\VMh + C\\B\\UW\h- 
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Now we turn to the equation for B. Taking t-derivative of (I3.29p 3 . multiplying 
by B t and integrating over Q, then 

-4 I \B t \ 2 dx+ [ \WB t \ 2 dx 
(3.72) 2dtJ 1 



= - J{u t - V)£ • B t dx + J{B t - V)u ■ B t dx + J {B ■ V)u t ■ B t dx 

Here Poincare inequality gives that 

- J \u t - V)B ■ B t dx 
(3.73) <|KIU4||VS|| £ 4p t || ia 



1 



<ollVu f ||! 2 + c||V5||^||5 t ||| 2 . 

o 

Gagliardo-Nirenberg inequality gives that 

J {B t • V)u • B t dx 
(3-74) <\\B t \\lA\^u\\ L 2 

<\\\Bt\\m+C\\Vu\\ 2 L 4Bt\\ 2 L 2. 

o 

And Holder's inequality gives that 

J (B ■ V)u t ■ B t dx 



(3.75) 



<7Tl|Vw t ||| 2 +C||S||2 || B .|2 



Collecting all the estimates f )3.63p -f l3T75]) and taking Proposition 13.21 I3.3j 13.41 
into account, we get that 



1 d 
2~dl 

(3.76) 



J \y/put\ 2 dx+^— J \B t \ 2 dx + ^ J \Vu t \ 2 dx + ^ J \VB t \ 2 dx 



<C(1 + \\u\\U + II^Hloo + + \\B\\ 2 H2 )(\\^u t \\ 2 L 2 + \\B t \\ 2 L2 ) 

+ C||Vu|| L 2||m||^2 + C||m||loo||Vu||^2, 

which together with Gronwall's inequality completes the proof of Proposition 13.51 

□ 

As a corollary, we can bound ||w|| L 2 W 2,4, which will play an important role in 
the estimates for p. 

Proposition 3.6. Under the assumptions of Theorem \1.3[ it holds that 
(3.77) sup {IMIz^cvtw 2 . 4 )} < oo. 

0<T<T* 
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Proof. It follows from Lemma 11.51 that 
(3.78) 

\\u\\iy2A 

< C\\u\\ m + C\\fm t \\ L 4 + C\\(pu ■ V)u\\ L 4 + C\\(B ■ V)B\\ L 4 

< C\\u\\ H i + C\\Vu t \\ L 2 +C\\u\\ L «>\\Vu\\ L * + C\\B\\ L «,\\VB\\ L i 

< c\\u\\ H i + c\\Vu t \\ L 2 + c\\u\\ L ~\\ VuW^WuW 1 ^ + c\\B\\ Loa livsii^iisn^ 2 , 

which finishes the proof of (13.771) . 



□ 



Furthermore, we have the following proposition. 
Proposition 3.7. Under the assumptions of Theorem \1.3\ it holds that 
(3.79) sup {||tt||#2 + ||5||_ff2} < oo. 

0<T<T* 

Proof. If the inequality (I3.48P is reconsidered, then the proof is done. 



□ 



Step V Estimates for ||V/>||Loo( 0>T ._ ff i) and \\(u, B)\\ L 2^ T . H 3y 
Proposition 3.8. Under the assumptions of Theorem \1.3[ it holds that 

(3.80) SUp (|MU~( ,T;if3)+ / (\H 2 H 3 + \\Bf H3 ) dt\ <oo. 
0<T<T* I JO ) 

Proof. Taking the Xj (j = 1, 2)-derivative of (13.291^ . 

(3.81) (p x .) t + u ■ Vp Xj = -u Xj ■ Vp. 

Multiplying the new equation by p Xj , integrating over Q, and summing up, then 
we obtain 

(3.82) — j \S7p\ 2 dx<C J \Wu\\Wp\ 2 dx < C\\Vu\\ L ^\\V p\\ 2 L 2. 

Similarly, we have the following higher order estimate for p, 
d 

(3.83) dt 



V 2 p\ 2 dx<C J (|Vw||V 2 p| 2 + |V 2 w||Vp||V 2 p|)dx 

< C||Vn|| L oc||V 2 p||i 2 + ||V 2 M |U4||Vp|| L 4||V 2 p|U 2 . 



Making use of Sobolev embedding inequality and Gronwall's inequality, we get 
that 

(3.84) HVpCOH^i < C||Vpo||iiexp Qf C\\ Vu{t)\\ w i,4d?J < oo. 

It follows from Lemma [1.51 that 

IMU 3 < C (\\u\\ H i + Hp^H// 1 + \\pu ■ Vw||#i + \\B ■ VB\\ H i) 

(3.85) < C(\\u\\ H i + ||Vp||La|K|| L a + \\u t \\ H i + || Vp|| L 2 \\u \\ L ~ \\ Vu\\ L 2 ) 

+ C (\\Vu\\ 2 L 2 + ||u|U=o||Vu||fli + \\B\\ 2 H1 + \\B\\ L oa\\VB\\ m ) 

which implies that sup 0<r<T , ||m||l 2 (o,t ; _h" 3 ) < oo. Similar proof leads to the same 
conclusion for B. This completes the proof of Proposition 13.81 □ 
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Combining all the estimates in Proposition 13. 2[ 13.51 and 13. 8[ we prove that 
( I3.28P holds and complete the whole proof of Theorem 11.31 
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